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Abstract 

The theory of degenerate parabohc equations of the forms 

ut = (^(mj;))^. and Vt = ($(f 

is used to analyze the process of contour enhancement in image processing, based on the 
evolution model of Sethian and Malladi. The problem is studied in the framework of 
nonlinear diffusion equations. It turns out that the standard initial-value problem solved 
in this theory does not fit the present application since it it does not produce image 
concentration. Due to the degenerate character of the diffusivity at high gradient values, 
a new free boundary problem with singular boundary data can be introduced, and it can 
be solved by means of a non-trivial problem transformation. The asymptotic convergence 
to a sharp contour is established and rates calculated. 
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1 Introduction. A model for contour enhancement 



This paper is devoted to study the behaviour of the solutions of some degenerate parabohc 
equations with moving boundaries which appear in describing the technique of contour en- 
hancement in image processing. Indeed, computer vision has become in recent decades a 
mathematical discipline which relies on the differential-geometric approach. More specif- 
ically, an appropriate technique of image processing consists of formulating a partial 
differential equation of evolution type for the image intensity, u{x,y). This function, also 
called the grey level, takes values in the interval < m < 1 and is defined on a two- 
dimensional image domain, Q. The usual evolution model leads to a nonlinear equation 
of parabolic type, possibly degenerate or singular. The nonlinearity is created by the law 
relating the image intensity flux to the image intensity. 

It has been observed by Perona and Malik |PMj . 1990, that for a suitable choice of 
this fundamental nonlinearity there appears an effect of enhancement of image edges that 
has a strong interest in the application to processing, denoising and recognition of images. 
The Perona-Malik Anisotropic Diffusion model has had a deep influence in the field, being 
at the source of many later developments. The model proposed by Malladi and Sethian 
|MSj leads (after proper scaling) to the following equation for the image intensity: 



(1.1) Ut = (i + \Du\'^y/^ K 

where Du denotes the spatial gradient of u, and k denotes the curvature of the surface 
z = u{x,y). The equation represents movement by curvature (curvature flow) and can be 
written as 

(1 + UI)U^^ - 2U^UyU^y + (1 + ul)Uyy 

l + ul + ul 

We consider here the more general flow given by equation ut = {1 + \Du\'^Y^^'p^/'^k, where 
p is a constant parameter. In other words, we study the equations 

(1 + 

with parameter a = —p/2. Along with the former case p = a = 0, the case p = —2, 
a = 1 has also attracted the attention of researchers (Beltrami flow, cf. Sochen et al. 

jHElS]). 

The asymptotic treatment of these models done in the papers |MS| ISKMSj shows 
the enhancement of the intensity contrasts by formation of regions of large intensity 
gradients, i.e., the normal component of the image intensity gradient becomes quite large. 
This phenomenon was analyzed in [BOlj . where a further simplification of the model 
was proposed in order to focus on the boundary layer where large gradients concentrate. 
Arguing locally around a sharp gradient point and choosing the the direction 
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normal to the boundary layer or front, we may disregard the effect of y derivatives with 
respect to the x derivatives in ()1.3p . In this way we get the reduced equation, which is 
just the one-dimensional version of ()1.3p 

(1.4) ut 



1 + uiy- 



where we have neglected Uy, Uyy. Different dimensional constants appear in the model, 
but they have been scaled to unity here without loss of generality. The mathematical 
problem consists in solving this equation with suitable boundary data, namely, m = on 
the left-hand side of the contour and u = 1 on the right-hand side (be that a finite or an 
infinite distance), and initial conditions 

m(x,0) = uo{x), 

satisfying < < 1 and Mq > in an interval / = (a, b) and constant values otherwise, 
zero to the left, 1 to the right. As was pointed out in |B01j . the phenomenon of gradient 
enhancement takes place in this model (in a proper setting) for all a > 0: the spatial 
gradient of the solutions, u^, increases with time, and its support shrinks. Indeed, we show 
below that such a behaviour can be observed in the larger exponent range a > —1/2. The 
conditions on uq can be relaxed, but then problems arise. Thus, the less stringent size 
restriction < mq < 1 or the lack of monotonicity create interesting alternatives, that will 
also be briefly discussed. 

Let us remark that this is not the only model which uses nonlinear diffusion equations 
for image processing and enhancement. We refer for alternatives to the works of Alvarez, 
Lions and Morel |ALMj and Caselles et al. |CMSj , where further references are found. Let 
us point out that in jALMj a fundamentally different assumption is made, namely that 
the flux is perpendicular to the gradient, and the basic equation is substantially different 
from the class ()1.3|1 treated here. On the other hand, we mention that the mathematical 
difficulties of the original model by Perona and Malik have been further investigated by 
several authors, like |Kchj and [WBJ. The application of Nonlinear Diffusion to Image 
Processing is a every active concern with many issues being discussed in the literature, 
cf. e.g. [KHIIW^ . 



2 Asymptotic self-similarity and enhancement 

Evolution equations like ()1.4j) and many other variants have been studied and are known 
in the literature under the general name of nonlinear parabolic equations of diffusion type, 
or nonlinear diffusion equations for short. They are typically used in describing processes 
of mass diffusion or thermal propagation. A quite general one- dimensional form popular 
among PDE experts is 

Ut = {a{x,t,u,Ux))x 
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with suitable conditions to make it parabolic, at least in a formal sense, like dajdu^ > 0. 
However, many practical applications (as in the present case) involve functions a for which 
da/dux > 0, but the values zero or infinity can also be taken, and then the equations 
are known as degenerate parabolic or singular parabolic resp., cf. jDB| IKaj . In any case, 
a general feature of this wide class of equations is the diffusive character, which roughly 
means the spreading of the level sets of the solutions as time advances. This property goes 
squarely against the desired enhancement, therefore an extra mechanism must be present 
if enhancement is to occur. We recall that in the Perona-Malik model this mechanism 
was negative diffusion. 

In the range of exponents a > 0, equation ()1.4|) falls into the class of degenerate 
parabolic equations with degeneracy at = oo, precisely the limit value which is of con- 
cern in contour enhancement. An investigation of the phenomenon of gradient enhance- 
ment in equation (jl.4|) is performed in the paper |B01j by relating it to the convergence 
toward self-similar asymptotics of an approximate model. It goes as follows: the observed 
evolution of the solutions towards a configuration with large gradients makes it plausible 
to further simplify the expression 1 + u^. in ()1.4|) into m^, so that the relevant reduced 
equation becomes 

(2.1) Mi = m;2{i+-)„^^. 

It is further observed that self-similar solutions for this equation with end-levels m = 
and u = 1 exist for all a > and exhibit the similarity form 

(2.2) m(x, t) = F(0, i={x- xo) it + to)'/'" , 

where Xq and are parameters to be fixed, and the profile F{^) is an increasing function 
joining the levels F = at a finite distance ^ = — c < to the level F = 1 at ^ = c. 
At these levels, which are taken at a finite distance ^ = ±c, the gradients are infinite 
(actually, the analysis in [BOlj deals with decreasing profiles joining u = 1 to m = but 
these two problems are obviously equivalent after a mirror symmetry. We have chosen the 
increasing option to avoid chasing around many minus signs). The scaling implies that 

u,{x, t) = {t + to)'/'"F'(0, i={x- Xo) {t + to)'/'" , 
which shows that the solution is concentrated in an increasingly narrower strip 

(2.3) ^ = {(x, t) : |x - XqI < c (t + to)"'/'"} 

with gradients that diverge like t^/^" as t — oo. Let us remark that the asymptotic 
divergence of the gradients is known to be exponential for a = 0, but the form is not self- 
similar of the same type, as we will see below. An important feature to be remembered of 
these solutions, along with the divergence as t — oo, is the infinite gradient condition at 
the end-points of the domain of definition. On the other hand, these self-similar solutions 
represent the intermediate asymptotics of the problem with more general data. This is 
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demonstrated numerically in |B01j when infinite flux data are imposed at moving end- 
points located at flnite distance. 

With these preliminaries we are ready to attack the construction and analysis of solu- 
tions with steep fronts for a general class of equations which includes ()2.1|) and ()1.4|) and 
is natural for our application, and with the general data mentioned in the Introduction. 
We use the theory of nonlinear diffusion equations which has been strongly developed in 
the last decades. It turns out that the standard initial value problem solved in this theory 
does not suit our model, since it leads to dispersion instead of concentration of the image. 
A new free boundary problem with singular boundary data is then introduced and solved 
by means of a number of non-trivial problem transformations available for nonlinear diffu- 
sion problems, a subject that is developed in detail in |V03cj . The concentration property 
of these solutions is possible thanks to the degenerate character of the diffusivity for high 
gradient values. 

Our theory below covers existence and uniqueness of solutions of suitable problems, 
existence and behaviour of the bounding interfaces, and large-time behaviour. In particu- 
lar, we obtain the rates of convergence of the interfaces, an important question for image 
enhancement. 



As we have just said, we are going to construct a mathematical theory for the one- 
dimensional evolution problem with initial and boundary conditions of the type mentioned 
above, but replacing the special class of equations ()1.4|1 by the wider class of nonlinear 
diffusion equations of the form 



The nonlinear function $ expresses the dependence of the image flux on the gradient of 
the image intensity. It can be called the flux law of our process, or also the constitutive 
function. In this general setting, $ can be any increasing real function deflned in a suitable 
interval of gradients. We will assume for convenience that 

(Hi) is C°° smooth and strictly increasing in the intervals < s < oo and 

— oo < s < 0, allowing for a quite arbitrary behaviour as s ^ or \s\ ^ oo. 

In our opinion, this level of generality instead of obscuring the problem makes more 
apparent the relation between nonlinear diffusion and contour enhancement. Slightly 
more general assumptions can be made on $, but they are unessential for our present 
purposes, cf. Section [TT] 

We want to characterize the class of those $ for which we can construct solutions 
with gradient-enhancement. We also want to describe the rates at which the formation of 
steep proflles takes place, thus justifying and extending the results of [BOlj to the actual 
equations like ()1.4|) . and to general data. The mathematical problem we pose consists in 



3 Nonlinear diffusion equations 
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solving the nonlinear parabolic equation (j3.ip with initial conditions 
(3.2) u{x, 0) = Uq{x) 

satisfying < mq < 1 and other suitable conditions (see below), and boundary data u = 
on the left-hand side of the contour and -u = 1 on the right-hand side. Regarding the 
initial data, we remark that we are mainly interested in monotone solutions, i.e., > 0, 
for equation ()1.4|) . and this condition will naturally follow from a similar monotonicity 
condition on the initial data. As for the boundary conditions, it turns out that, depending 
on the form of $, the boundary can be chosen to be located either at infinity or at a finite 
distance. This latter case will be the one of interest for us, and then the problem must 
be properly posed as a free- boundary problem. 

Let us recall the simplest examples of $ and some of the difficulties we will encounter. 
Indeed, a quite important and simple example is the power function, that we write as 
$(s) = (l/m)s™, defined for s > 0, so that it agrees with 1)2.11) with m = — 1 — 2a since 
$'(s) = s"^~^. There is in principle no reason to restrict the generality of the exponent m 
in the mathematical treatment to follow, and this will lead to quite different behaviour 
types inside this family. The case m = is included in the form $(s) = log(s), i.e., 
$'(s) = Note that in the cases with m < the function $ is negative, but the 

important quantity for the parabolic character of the equation, i.e., ^'{ux), is always 
positive. Finally, we note that equation ()1.4p corresponds to $'(s) = (1 + s'^)"^^'^"^; it 
degenerates as s — > oo in the prescribed range a > 0, even if a > — 1, but the equation is 
is perfectly parabolic in regions of bounded Ux- 



3.1 Non- monotone solutions 

The above statements are made on the assumption that the solutions are monotone, 
Ux > 0, which is not unjustified in our problem setting but restricts the mathematical 
generality. When considering non-monotone solutions it is customary in the nonlinear 
diffusion literature to extend the power nonlinearity to arguments s = m^,. < in the 
simplest symmetric way: 

(3.3) <l>(s) = — Isp-^s, 

m 

which gives $'(s) = |s|™~^, always nonnegative. This definition poses no problem when 
m > and solutions with changing sign exist corresponding to initial data with the same 
property. However, we are interested in exponents m < (so-called very fast diffusion 
in the literature) where $' is singular at s = 0, and the whole function $ is no more 
monotone. The difficulty has been studied in |RV02j leading to the consequence that 



solutions with changing sign do not exist, even in the weak sense when m < 0. This 
is a consequence of the singularity of the equation at the level Ux = and does not 
affect equation (|1.4p . Due to this obstruction, we will concentrate here on problems with 
monotone solutions. 



6 



4 Second formulation as a nonlinear diffusion equation 



If we formally differentiate equation with respect to x and put f = Mj;, we obtain 
the equation satisfied by the image intensity gradient: 



which is usually called the nonlinear filtration equation, NLFE, and is the most standard 
class of nonlinear diffusion equations studied in the literature. We can call it in this 
context the "differentiated equation" . Inversely, we can recover the u formulation from a 
solution v{x,t) for ()4.H) by means of the rule 



integrated along any curve T in the domain of definition of v which joins a fixed point 
(xo,to) to the generic point {x,y), cf. |RV95[ IV08cj . The constant c is the value of u at 
{xo,to) to be chosen at will in principle. The calculation is justified for smooth solutions 
V and smooth $, but holds in a much wider context. 

Monotone solutions u for equation p.4j) translate into nonnegative solutions v for 
equation ()4.1|1 . and conversely. In this context, u is usually viewed as the mass function 
for V, since when we take F to be a segment of the line t = t^, formula ()4.2p becomes 



The phenomenon of gradient enhancement can be then translated into usual diffusion 
parlance as mass concentration. We will keep in the sequel the denomination intensity or 
image intensity for the solution u of equation (jH.lj) . and we will view the solution v = Ux 
of 1)4.11) as the intensity gradient. In that context ^{ux) = $(f ) is the intensity flux. 
Finally, = ) is the diffusivity. 



We now try to apply the more standard theory of nonlinear diffusion equation to solve 
the evolution problem motivated by the application in the Introduction. We show next 
that this implementation can be performed in a rather standard way for a wide class of 
$'s, but such a process gives an evolution with no concentration effect, hence useless for 
our purposes in image processing. The reader already familiar with diffusion theory may 
check the contents of Theorem 15.11 and skip the rest of this section. 



(4.1) 



(4.2) 





5 Solutions of Type I 
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5.1 Simplest setting 



Generalizing the well-known properties of the heat equation (the choice = s), we 
consider first the case where $ is a function and does not vanish. We pose the 
Cauchy problem for equation ()3.1|) in the whole line x G M for t > with bounded initial 
data < Mo < 1, and get a unique smooth solution u{x, t) defined in Q = ffi x (0, oo) and 
such that < u < 1, M is smooth for all a; G M and t > r > 0. We can also work with 
equation (jUIJ: then, if vq = uo,x is locally integrable (a Radon measure will also do) and 
such that 



the solution v{x,t) is smooth in Q and satisfies the same type of integrability condition, 
J^v{x,t) dx = 1 . For solutions f > this is called conservation of mass. Using formula 
fl4.2|l we obtain a smooth solution u{x,t) in Q such that the initial condition 



hold locally uniformly in time. On the other hand, f o > implies that Ux = v > 0, which 
by the Strong Maximum Principle implies > 0. It follows that < u{x,t) < 1 in Q. 
Note that, due to the smoothness condition on $, these solutions are smooth in {x,t) for 
all X G M and t > 0, even if the initial data are not. 

5.2 Dispersion of solutions with time 

At first glance, the type of solution we have constructed seems to solve our image problem. 
However, it lacks a basic ingredient, i.e., the eventual concentration of intensity gradients. 
On the contrary, if we consider the heat equation Ut = u^x, the solutions, which can be 
expressed in terms of error functions, spread in time and its gradients go to zero, for 
instance, if Vq = Uq^x is integrable, then v{x,t) goes to zero as t — oo at a rate of the 
order of and takes on a Gaussian space shape, 

v{x,t)f/'^ ~ exp(-xV4t). 

And a similar result (with possibly different rates) applies to more general data uq having 
definite limits at ±oo, and to all the functions $ of the above class. 

5.3 Extension of Type I solutions to other $ 

Even if this is not the class of solutions we are looking for, we pursue a bit further the 
analysis, since it is the standard class found in studies of nonlinear diffusion or thermal 





holds, as well as the end conditions 



lim u{x,t) = 0, lim u{x,t) = 1, 



— oo X— >oo 



8 



propagation. The study will also serve for comparison with the "correct" solutions of 
Type II. However, the reader may choose to skip this the rest of the section and proceed 
with Sectional 

The class of constitutive functions $ for which there exists a class solutions < m < 1 
with a dispersive character can be extended to include a quite general choice of $ on the 
condition of allowing for a suitably generalized concept of solution. Thus, it is well-known 



following Benilan, Crandall and others |Be72[ IBC'j that we may take as $ any continuous 
nondecreasing function and the Cauchy problem is then well-posed in the class of so-called 
mild solutions with L^(M) data. Actually, $ can be allowed to be discontinuous but this 
will be of no particular interest here; on the contrary, we will keep the assumption of 
smoothness and strict monotonicity for < s < oo for simplicity of presentation (and 
since it is satisfied in the application we are dealing with). 

Let us comment on the main properties of the solutions for this class of equations. For 
the case of power nonlinearities $(s) = s'^/m, s > 0, mentioned above, rather complete 
details are known, |Ar| IKa| IV92j . We can consider that m is a real parameter, positive in 
principle but not necessarily as we will see. The definition is extended to solutions with 
negative values as in formula (jH.Hj) . Thus, for m > equation vt = {v"^)xx (or better, 
Vt = {\v\''^~^v)xx) generates a positive semigroup of contractions in the space L^(M). In 
other words, for every Vq G L-'^(]R) there exists a unique function v G C([0, oo : L^(M)), 
f > with v"^ G L\^^{Q) such that the equation is satisfied in the sense of distributions 
in Q = M X (0, oo), the initial data are taken in the L-'^(]R)-sense, and the map vq ^ f (■, t) 
is an L^-contraction. Moreover, the total mass is conserved 

v{x,t)dx= / VQ{x)dx. 
Jo 

In our application we still have to impose the extra condition of total mass 1, and then 
we recover the image intensity by means of the formula 

u{x,t) = I v{x\t) dx' , 



and the intensity level goes from m = at minus infinity to u = 1 at infinity. The 
asymptotic behaviour of these solutions has been carefully calculated in the literature, cf. 
e.g. |Xn2llFKHVn8aj : solutions t;(-, t) with finite mass go zero uniformly as 

(5.1) v{x,t)^r^, 7 



m + 1 

which in the notation of the Introduction means 7 = —l/2a. Actually, the asymptotic 
rate comes from comparison with the source-type self-similar solutions |Bbkj which take 
the form 

v{x,t) =t-^G{xt~^) 

for a certain symmetric profile function G G Cfe(M) such that G(0 ^ as ^ — 00. This 
formula suggests that the behaviour will be maintained as long as 7 > 0, hence, as long 
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as m > —1, and this turns out to be true (for the study when — 1 < m < cf. |ERVj ). 
Indeed, this exponent range is optimal since there are no solutions with finite integral for 
m < —1. 

Let us look a bit closer at the kind of initial data that we want to consider in the 
outmost generality. We recall that we want Uq to be increasing and bounded between 
and 1. Hence, Vq = u'q has to be defined and locally integrable, or at most be a measure, 
in an interval (a, b) with possible divergence at the end-points, but anyway with finite 
integral. Existence of solutions (with vq extended all of M with value otherwise) offers 
no problem in all the range m > — 1 and the solutions are bounded for alH > 0. But large 
values will be the origin of the new class of solutions to discuss in the following sections. 

By integration we obtain a solution of the problem 

( ut = ^{u^)^ in g = Rx(0,oo) 

(P/) < u{x,0) = uo{x) for xeR 
[ ulx,t) > 0, 

where Uq is any nondecreasing continuous real function joining the levels Mq = at x = 
—oo to uq = 1 at a; = oo. This is what we term as Solutions of Type I. Note that 
integration of the self-similar profile gives for u the form 

u{x,t) = F{xr^) 

where F is a primitive of G with F(— oo) = 0, -F(oo) = M that we want to normalize to 
M = 1. 

Theorem 5.1 Let us consider the initial-value for equation (j3.H) posed in Q = 'Rx {0, oo) 
with power function $(s) = s'^/m. If m > —1 then for every nondecreasing uq with 
Mo(— oo) = 0, Mo(oo) = 1, there exists a unique continuous weak solution in the sense 
u{x,t) > such that u{-,t) jumps from to 1 as x ranges over the line x G M. If m < 
the last condition is essential to ensure uniqueness. This class of solutions has bounded 
gradients for strictly positive times (t > r > 0), and spreads in space as time advances 
and the asymptotic decay formula 1)5.11) holds for v = u^. 

On the contrary, if m < —1 solutions for this initial-value problem with bounded data 
do not exist. 

The results of the power case refiected in this theorem can be generalized to equations 
with more general Only the behaviour of $(s) at s = and s = oo will determine the 
different behaviour types. In order to tackle the first, we assume that the initial data (and 
hence the solutions) v are bounded. Then it is known that the condition for existence 
with finite mass is 

(5.2) [ {s)s ds <oo. 

Jo 

The question of large arguments is similar to the power case. These issues will be discussed 
at length in )Vbkj . For a study of these questions in the class of self-similar solutions cf . 
|FV|. 
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5.4 Existence of sharp interfaces 

There is an interesting subclass of equations (jH.lj) . or (jHH), where sharp interfaces appear. 
Let us look first at power nonlinearities. Indeed, for exponent m > 1 solutions of the 
Cauchy problem for equation ()3.1|) with initial data having compact support will keep this 
property for all times, while infinite propagation occurs whenever m < 1 cf . |BV| IKaj . 
In the first case, given an integrable function uq > with integral 1 the solution of the 
problem can be seen as the a classical solution equation (jH.lll in a domain 

n = {{x,t) : -l{t) <x < r{t)} , 

with initial conditions u{x, 0) = uo{x) and boundary conditions 

( u{x,t) = 0, ^(u-^) = for x = l{t), 
t) = l, <l>(-u^) = for x = r{t). 

The lines x = l{t) and x = r{t) are called interfaces or moving boundaries and are 
completely determined by the over-specified conditions. They are known to be smooth 
(analytic) functions of t |Anll lAVj . and diverge as t cxo like 0{P), thus giving a 
quantitative estimate of the dispersion effect. We remark in passing that the presence of 
interfaces means also that the equation is not uniformly parabolic at those points, and 
consequently the solutions have limited regularity. 

On the other hand, when m < 1 the same Cauchy problem leads to positive solutions 
with l{t) = — oo and r{t) = oo. The property of null flux is equivalent to imposing Ux = 
at ±cxD, a quite natural condition in view of the values m = 0, 1 at ±oo. This condition 
is automatic for m > 0. However, for m < we can have new solutions with decreasing 
total mass, i.e., such that 

— / v(x, t) dx < 0, 
dtj ^ ' ' 

and they can be even forced to vanish identically in finite time by controlling the outgoing 
flux at X = ±cxD. We refer to |RV93|rRV95j for a detailed analysis. In the integrated version 
they would lead to solutions with a restricted grey range. For the more general class of 
$ mentioned above the condition to have finite interfaces is 

(5.3) / ^^ds<oo. 

Jo s 

The fact that this condition is necessary and sufficient can easily seen on the family of 
traveling waves. The property was first pointed out in |()K(yj . 

6 Solutions of Type II. Conjugate formulations 

We now address the main question of this paper, the construction of solutions with large 
gradients, appropriate for the contour enhancement problem. 



11 



6.1 Basic Free Boundary Problem. It is formulated as follows: 
Given an increasing function Uq{x) defined in an interval (a, b) with end values u{a+) = 0, 
u{h—) = 1, to find a continuous function u{x,t) and continuous curves x = l{t) and 
X = r{t) such, that 

(i) /(O) = a, r(0) = h, and l{t) < r{t) for some time interval t e (0, T). 

(ii) u solves the following problem in Q — {{x, t) : < t < T, l(t) < x < r{t)}: 



Such a triple {u,l,r) is called a solution of Type 11. The regularity required from u as 
a solution of the problem will depend on the generality of the data. At least u will be 
continuous in the closure of Q. Furthermore, to avoid unnecessary generality we will ask u 
to be smooth in the interior of Q. Finally, the requirement of monotonicity is not intrinsic 
from the mathematical point of view, but it suits the application and allows for the use 
of the powerful conjugate formulations. Here is the existence result that we are going to 
prove: 

Theorem 6.1 Let ^ be a flux function that is defined, smooth and with > for all 

s > 0. Assume moreover that $(00) is finite. Then for every increasing function Uq{x) 
defined in an interval [a, b] with u{a) = 0, u{b) = 1 and Mq > c > there exists a unique 
continuous function u{x,t) which is defined in a set fl as above, is smooth and strictly 
monotone in x for < u < 1, and there exist continuous curves l{t) and r{t), such that 
the triple {u, I, r) solves problem (Pjj) in Qt. Besides, Ux > c> whenever < u < 1. 

We will also show in our construction that the boundary curves are monotone: l{t) is 
nondecreasing, r{t) nonincreasing. On the other hand, there is the problem of determining 
whether T if finite or infinite. This depends on $ as we will see below. 

6.2 Conjugate formulation 

When dealing with smooth monotone solutions m^: > we can invert the variables x and 
u and write x = X{u, t). Then u^-x^ = 1, and after some computations we get the partial 
differential equation satisfied by function of u and t: 



(Pii) 



I 




u{x, 0) = uo{x) for a < X <b 

u{l{t),t) = 0, Ux{l{t),t) = +00 for 0<t<T, 
u{r{t),t) = 1, Ux{r{t),t) = +00 for 0<t<T. 



(6.1) 



Xt = 



where ^ is the conjugate flux function (conjugate to $), defined for s > as 



(6.2) 



$(l/s) . 
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Differentiation of equation with respect to u gives rise to tlie differentiated conjugate 
equation for w = dX/du clS 8b function of u and t: 



(6.3) wt = {^{w))uu- 

We complete tlie list of related equations with the direct differentiated equation for v = 
du/dx, already seen: 

(6.4) vt = (Hv)),, , 

and then v = 1/w. It is important to point out that these relations are equivalent to the 
well-known Backlund transform |BKj between the main variables v and w of the second 
and fourth formulations. Indeed, we have 

1 r 

v(x,t) = — r, u(x,t)=c+ (v dx + ^(v)x dt) . 

w{u,t) Jr 

The reader is referred to |RV95| IV03c| IV03dj for other applications of this technique. 
6.3 Posing and solving the conjugate problems 

We now show how to use the conjugate formulations to solve the original problem. We 
assume that $ a flux function defined for all s > and such that $(oo) is finite, say, 
$(00) = 0. Then \E'(s) is continuous at s = and ^^(0) = 0. This is the class of flux 
functions for which the conjugate problem looks simpler. Since we have assumed that $ 
is smooth, so is in its domain of definition. 

(i) Since we are interested in solving the conjugate problem as an auxiliary step for 
Problem (Pjj) we will relate the initial values for the function w{u,t) to uq as follows. 
Assuming that uq is continuous and strictly monotone in the interval I = {a < x < b}, 
with uo{a) = 0, uo{b) = 1, and smooth inside I with duo/dx bounded below away 
from zero, we define the inverse function x = h{u) = Uq^{u) : [0, 1] [a, b]. 

1 

(6.5) Wo{u) 



which is defined for < u < 1 and is positive, bounded and smooth inside, i.e., for 
u G (0,1). 

(ii) We then solve the conjugate problem 

{Wt = '^iw)uu for < M < 1, it > 

w{u,0) = woiu) for < M < 1 
w{u,t) = for M = 0, 1. 

As initial data we choose a nonnegative, bounded function wq. Under these conditions 
Problem (P^) has a unique solution by virtue of well-known nonlinear parabolic theories 
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described in the previous section, cf. |Be72| \nC\ IBGj : but note that now we are deahng 
with the homogeneous Dirichlet problem. The solution can be obtained as limit of the 
solutions Ws{u, t) > e of the nondegenerate problems with initial data tyo,e(u) = WQ{u)+e, 
6 > 0. In the monotone limit we get 

\imws{u,t) = w{x,t) 

which is nonnegative, continuous and bounded. Under the additional assumption that wq 
is locally bounded away from zero, it is easily proved that the solution w{u, t) is positive, 
hence classical, in a strip 

St = {{u,t): < M < 1, < t < T}. 

(iii) We also need to know something about the large time behaviour of the solutions to 
this problem. In the full generality (\E' continuous at s = 0), it can be proved rather easily 
that w{x, t) goes to zero in uniform norm as t — > oo. However, the rate depends on \Ef, as 
we will see in detail in the next sections. Indeed, depending on ^ it may happen that the 
solution vanishes identically after a finite time T > (so-called finite-time extinction). 
For power nonlinearities "^'{s) = this happens iff < g < 1, |BHj . On the other 
hand, it is well-known that for q > 1 the decay rate is 0{t^^^^'^^^^), while for g = 1 it is 
exponential. 

(iv) Next, we pass to the integrated version using the formula 

(6.6) z{u,t) = J w du + '^{w)udt, 

where F is any piece-wise smooth curve in {u,t) space starting from a fixed point, say 
u = l/2,t = and arriving at a generic point {u,t). In this way we obtain a solution of 
the integrated equation Zt = ($(2;„))u, much as we did in the case of the original pair of 
formulations. Note that 

(6.7) z{l, t) - z(0, t) = j wdu-^Q 
as t ^ T. 

(v) Let us examine the curves Z{^u){^) = z{u,t) for fixed u G [0,1]. It is clear from the 
smoothness of the solutions that these curves are C°° smooth for every < m < 1. We 
are interested in the limit curves 

(6.8) l{t) = \imz{u,t), r{t) = \im z{u,t). 

The limit is well-defined for every < t < T by monotonicity of 2 as a function of u. 
These curves will show up in the next subsection as the interfaces of the original problem. 
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Lemma 6.2 The curves z = l{t) and z = r{t) are continuous and monotone for all 
< t < T with a < l{t) < r{t) < b (r{t) is nonincreasing and l{t) nonincreasing) . As 
t^T we have r{t) - l{t) 0. 

Proof. For < s < t < T and < Mi < ^2 < 1 we have 

z{ui,t) — z{ui, s) = / w{u,t)du+ / w{u, s) du + {z{u2,t) — z{u2, s)). 

J Ul J Ul 

Fix now s > and let t be a bit larger than s. Since w is bounded uniformly for t > s 
the two integrals are uniformly small as long as ui — U2 is small. We will fix now M2 ~ 
and let ui G (0, U2) go to 0. It is clear then that z{ui,t) — z{u\^ s) is uniformly small and 
goes to as t ^ s. Hence, lif) is continuous at t = s. The argument for rif) as m — > 1 is 
the same. 

For classical solutions the monotonicity of the limit curves is a consequence of the 
estimate 

dz d 

m = 

which must be positive at m = and negative at m = 1 because > and we have zero 
boundary conditions. For the general case we use the dependence of the solutions on $ 
as demonstrated in |BCj . Hence, l'(t) > 0, r'(t) < 0. We will see another proof below. 
Note that the rest of the curves Z(u)it) need not be monotone unless wq is a rearranged 
symmetrical function. 

As for the last statement, it follows immediately from 16.71 It will mean for the original 
problem that the solution concentrates into a vertical profile, thus showing the formation 
of the desired front. ■ 



6.4 Inversion. Solutions of Type II 

Thanks to the fact that dz/du = w > 0, we can invert the dependence between z and u 
in the previous construction to get a function u = u{z,t) that is easily shown to satisfy 
the equation 

Ut = 

Besides, m is a monotone function of z and takes the values u = and u = 1 respectively 
at the left and right endpoints of the domain of definition 

= {l{t) <z< r{t)}, l{t) = z{0, t), r{t) = z{l, t). 

where z{-,t) is the function defined in ()fi.fi|l . Therefore, u{z,t) is a candidate to solve our 
original problem if we identify the independent variable z with x — c, where c is uniquely 
determined by the relation mo(c) = 1/2. 
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In order to check that we have solved the original problem we still have to check some 
particulars. It is clear that v = is related to the original w by the formula 



v{x, t) 



1 



w{u, t) ' 



which simply states the derivative rule for the inverse function, and u in the second 
member is given by u{z,t), explained before. Here comes an important 

point: since w takes on zero boundary values, v{x,t) diverges at the endpoints of its 
domain of definition, Q. In other words, the solutions of the original problem u = u{x,t) 
enjoy the property of infinite gradients at the endpoints of the strip where they are defined. 
Since $(cx3) = this also means zero flux at these points, a reasonable requirement, which 
explains why this condition has to be imposed on $. 

As for the initial data, we have the mass formula 



so that X ranges over an interval [a,b] when u goes from to 1, i.e., a = -2(0,0) + c, 
b = z{l, 0) + c. This rule is accompanied by the rule u = J^^^-j v{x, t) dx. 

6.5 Uniqueness 

Uniqueness of our class of monotone solutions works by translating any couple of solu- 
tions of Problem (Pjj) with the same initial data into the conjugate formulation. They 
continue to have the same initial data. Uniqueness of weak solutions is well-known for 
that equation. The Theorem is proved. ■ 

6.6 Front formation 

The asymptotic formation of a vertical front is a simple consequence of the fact that there 
exists the limit 



7 Uniqueness and comparison for Type II solutions 



The previous construction provides for existence and uniqueness of monotone solutions. 
We tackle next the property of comparison, which is stronger that uniqueness. 




1/2 




1 

2' 



lim r(t) = lim l(t) = Xoo ^ (a, b). 
The existence of the common limit follows from Lemma [6.21 



16 



Theorem 7.1 Comparison applies to the solutions of Theorem Va.W : letui{x,t), i = 1,2, 
be two solutions defined in a strip 5^ = K x (0,T) having initial data Uoi, where the 
solutions and the data have been extended by the the left of the definition domain, by 1 
to the right. If Uoi{x) > Uo2{x), then for a// x G M we have 



Proof. The main idea is to use a contradiction argument at touching points of ordered 
solutions. We proceed in several stages and need a definition, taken from the theory of 
viscosity solutions in |(]Vj . We say that two solutions are strictly separated at a time t if 
Ml > U2 at all points intermediate to the bounding interfaces, and these are also separated. 

(1) Given two solutions Ui and U2 with ordered data we must obtain approximations 
to which the contradiction argument can be applied. First, we separate the initial data of 
the solutions by displacing the first solution ui to the left (the problem is invariant under 
space displacements). Next, assume that we have modified Ui so that it has a finite slope 
on the left-hand free boundary, and U2 so that it has a finite slope on the hand-side free 
boundary. In this situation we argue on the first point and time {xq, to) where the graphs 
of both solutions touch and discover that: (i) it cannot happen with value u e (0, 1) by the 
Strong Maximum Principle; (ii) not on the left-hand side where Ui = U2 = 0, because at 
such a point we have the right derivative: U2^x = +oo while Ui^^ is finite, which contradicts 
the fact that we still must have Ui(-,to) > U2{-,to)] (iii) same argument applies to the 
right-hand side where Ui = U2 = 1. We conclude that such solutions Ui and U2 cannot 
touch, hence strict separation is preserved in time, and comparison ui > U2 holds in the 
strong sense, i.e., with strict separation at all times. 

(2) We must now prove that the original solutions can be approximated by solutions 
with finite derivatives on the lateral boundaries. This can be obtained by an easy mod- 
ification of the previous construction. Thus, if u is the constructed solution, it can be 
approximated at the level of the conjugate problem by putting the value w = e on the 
left-hand side boundary m = 0. It is rather standard monotonicity argument that We ^ w 
in the limit e — > 0. Undoing the transformation this means that Ue converges to the orig- 
inal constructed solution u. A similar construction applies to U2 by modifying W2 at the 
border u = 1. After displacement of Ue we get + 5, t) > M2,e(a;, t). Let now e — >■ 0. 
■ 

We derive next another proof of the monotonicity of the interfaces. 
Corollary 7.2 The interfaces are monotone in all cases. 

Proof. Given a solution u with initial interfaces a and 6, we may place a very steep solution 
u to the left, i.e., with -u = for x < a — 2e, u = 1 for x > a — e. If -u' is symmetric then 
the interfaces are monotone, hence the left interface of u lies to the right of a — 2e. Then 
u > u for all times, hence l{t) > a — 2e, and in the limit l{t) > a. The same argument 
applies taking origin of times at any y G (0,T), hence l{t) is monotonically increasing in 
(0,T). Likewise, r{t) is proved to be monotone decreasing. ■ 
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8 Self-similarity and asymptotics for power nonlinearities 



In the power case, where = (l/m)s™', the asymptotic condition "$(00) finite" means 
m < 0. Then = {l/q)s'^ with q = —m > 0. Reminder: we consider only monotone 
solutions of ()3.1|) and nonnegative solutions of ()4.H) . and g = 1 + 2a; in the notation of 
the first section. Therefore, a > —1/2 means g > 0. 

There is a unique solution of problem (P^) with initial data as in Section IHl The large 
time behaviour depends on the exponent. For q > 1 solutions are defined and nontrivial 
globally in time, while for < g < 1 they exist only in a finite interval < t < T and 
w{x, t) — > uniformly as t — > T. 

As for the exact behaviour, if a 7^ the self-similar behaviour reduces to writing 
U = <l>(^) with ^ = xt^/^" and integrating the ODE 

(8.1) $"=^^($')2+9 

2a 

to obtain in the monotone region 

(8.2) $'(0 = AiK'~ ^2)-V(2+2a)^ ^ ^ (2a/(^ ^ a))i/(2+2") . 

and a constant K > determined by the end values min $(0 = 0, max $(0 = 1. For the 
large-time behaviour of general solutions we have to make a study in different cases. We 
start with a = which was not included before. 

Case a = 0, hence q = 1. The conjugate equation is the linear heat equation. It is 
well-known that the solutions behave for large times like 

(8.3) w{u, t) ~ Ci W{u, t) = Ci e-^*/i(M), 

where A = vr^ is the first eigenvalue of the Laplacian in [0, 1], /i(x) = sin(7ru) is the first 
eigenf unction, and Ci > is a constant that depends on the initial data. Here and in the 
sequel we denote by capital letters the quantities corresponding to self-similar solutions. 
We get the intermediate step 

(8.4) x-c = z{u, t) ~ Z{u, t) = -C e"-^* cos(7ru) , 

with C = Ci/vr. This formula can be viewed as an implicit expression of m = u{x,t) 
defined in the space between the interfaces. Note that M = l/2fora; = c and all t > 0. 
Putting c = we get for the self-similar solution 

(8.5) U{x,t) =F{xe^''/C), 

where F' behaves like d~^^'^ near the end-values or interfaces, d being the distance to these 
points. Moreover, the interfaces 



(8.6) R{t) =CRi{t), L{t) = -CRi{t), with Ri{t) = e 
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These estimates become the first order approximation when we consider general solutions: 



(8.7) M(x,t) ~ F(xe-" VC), r(t), /(t) ~ c ± C i?i(t). 

The gradients are given by 

1 1 



18) u^(x,t) = — ^—e 



xt 



w C sin(7rM) 

They blow up exponentially as t ^ oo and like Oi^d"^/"^) at the interfaces, where d is 
distance to the interface. 



Case a > 0, hence g > 1. This is the case treated in |B01j by the direct self- similar 
method. We recover the behaviour of general solutions from the conjugate problem. We 
have a self-similar asymptotic expression in the form of separation of variables: 

(8.9) w{u, t) ~ W{u, t) = r^/(^-^)/,(n), 

where fq{u) > is the unique solution of the associated elliptic problem 

(/T + /^/ = in [0,1], /^=^, 

and / = at the endpoints. This problem was studied in |APj . In this case there is no 
free constant, i.e., the behaviour is universal. From it we get the intermediate asymptotic 
estimates 

(8.10) x-c = z{u, t) ~ t-^/(''-^)^q(u) . 

This formula defines implicitly u = u{x,t) in the space between the interfaces as 

(8.11) M~Fg((x-c)t^/(«-^)). 

Since f{u) behaves exactly like 0{u^^'^) near the end point u = 0, and 0((1 — uY^'^) near 
M = 1, it follows after integration and inversion that u = F behaves like 0(rf(^/(''+^)), 
where d is the distance to the interfaces. These are given by 

(8.12) /(t) = c-c^r^/^"-^), r(t) = c + c^r^/^^-i). 

The gradient is given by 

1 



(8.13) M^(x,t) -t^/^"-^) 



and blows up like Oit^^"^") as t — > cx) and like 0{d '>') at the interfaces, where 7 
l/(g+l) = l/(2a + 2). 
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Case — 1/2 < a < 0, hence < g < 1. This is a new case, not included in the 
modehzation of |MS| IBOlj . We get solutions from the conjugate problem which exist for 
a finite time T and behave as t T like the separation of variables formula: 

(8.14) w{u,t)^{T-tf^'-'^'>f^{u), 

where fq{x) > is the unique solution of the associated elliptic problem 

+ /i/ = in [0,1], /i ^ 



and / = at the endpoints, with no free constant. We get by integration 

(8.15) x-c = z{u, t)^{T- tf/^^-'i'^gqiu). 

This formula defines implicitly u = u{x,t) in the space between the converging interfaces 

(8.16) l{t) = c-C{T -ty/^^-''\ r{t) = c + C{T-ty/^^-'^\ 
The gradient is given by 

(8.17) «^(a;,t)^(r-t)-va-^)-i^, 

fqW 

and blows up like 0{{T — t)^/^") as t ^ T and like 0{d~'^) at the interfaces, 7 as before. 



9 Asymptotic estimates for more general ^ 

Assumptions on $ to behave like a power at infinity guarantee that the asymptotic be- 
haviour is as predicted in the power case. Thus, Bertsch and Peletier study in |BPj the 
asymptotic behaviour of equation (jfi.lj) written in the form 

P{v)t = v^^. 

In our application the space variable x becomes u and (3 = \1'^^. They assume that 
/?(0) = 0, (3'{s) > for s > and 

-l<£^<a 

for all small < s < Sq and some a G (0, 1). They also need the more stringent condition 
that the limit 

a(s) = lim ^ , , 
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exists for every s > 0. This limit is necessarily of the form a{s) = for some m G [0, 1). 
Those conditions cover in particular powers "^{s) = s'^ with g = 1/m > 1, as well as the 
exponential '^{s) = e~^/**, and any other function resembling such examples near s = 0. 

Under these restrictions, they establish an asymptotic separation-of-variables result 
that, when translated to our setting, means that the solution w{u,t) of the Cauchy prob- 
lem ()6.1|) behaves as t — oo in the separate variable form w ~ y{t)f{u), more precisely. 



(9.1) hm ^ = f{u), 

t^co y{t) 

where y{t) is the solution y' = —y^^'^, i.e. y{t) ~ ct~'^^^'^~^^ and / is the profile obtained 
as in the previous section. This means that the analysis of the previous section is justified 
for almost power-like flux functions if they resemble a power with exponent g > 1. The 
cases q = 1 and g < 1 should be justified in a similar manner. We note that a similar 
calculation for the Cauchy problem was done in |PVj . 

We point out that these restrictions on the nonlinearity are satisfied by the equation 
proposed in the image processing model, ()1.4p . hence the results apply to that model. 



10 Numerical experiments 

We have computed the solutions of the free boundary problem (Pn)) by solving the 
conjugate problem, which is a homogeneous Neumann problem, in the cases a = 0, a = 1 
and a = —1/2. The conjugate equations are p-Laplacian equations with p = 2,4 and 
3/2. The computations have been done with symmetrical data Xo{u) = —cos{nu) and 
asymmetrical Xo{u) = 2u^ — [l/2)u^ — 1. The last two examples stabilize in finite time 
(since a < 0). 




Figures 1 and 2. The case a = with symmetrical and asymmetrical data 
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Figures 3 and 4. The case a = 1 with symmetrical and asymmetrical data 
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Figures 5 and 6. The case a = — 1/2 with symmetrical and asymmetrical data 

11 Discussion and conclusions 

We have established the well-posedness of a free boundary problem that represents a 
one-dimensional version of the model for image contour enhancement. We have also 
established the asymptotic behaviour and its rates. The results apply to a large class 
of equations, which enlarges the scope of the results of |BOH IMSj . We have performed 
a classification of the solutions and their properties according to the properties of the 
constitutive function $. 

We have used as a framework the theory nonlinear parabolic equations of diffusive 
type, or nonlinear diffusion equations for short. These equations are typically used in 
describing processes of mass diffusion or thermal propagation (with or without additional 
effects, like convection or reaction). Here they appear in image processing. The types used 
here cover the model cases known as p-Laplacian equations and porous medium equations 
or their variants. 
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The analysis is technically performed by means of a series of remarkable transforma- 
tions that lead to conjugate problems that are easier to analyze. These transformations 
are related to the Backlund transform. 

There are a number of interesting consequences of the mathematical analysis that we 
sum up next. 

• Generality. The first one is the observation that solutions with gradient blowup 
can be obtained for the nonlinear diffusion equation with a large class of constitutive 
functions. Namely, $ can be any continuous nondecreasing function defined for s > 
sucli tliat $(cxd) is bounded. We have made for convenience the assumption of smoothness 
and strict monotonicity for < s < oo. Actually, the constitutive function can be more 
general, even a discontinuous maximal monotone graph as in |Be72| IBrzj . But such a 
generality is not practical here, though the classification and results go through after 
some heavy work. A general presentation of that generality will be the object of the work 



• Focusing in finite time. In the case of powers ~ s"* the condition means 
m < 0, and not only m < —1 as considered before. Moreover, the new range > m > — 1 
leads to a very fast evolution which arrives at a vertical front in finite time. This idea 
might have an applied interest. 

• Coexistence of types. We have also a so-called class of solutions of Type I that fall 
into the scope of the standard parabolic theory. They exist under the condition that the 
integral J ds converges at s = oo, which in the case of powers means m > — 1. 
This fact leads to an interesting observation: there is a range of co-existence of solutions 
of Type I and II for — 1 < m < 0, i.e., —1/2 < a < 0. This may be of interest from 
the theoretical point of view: two different problems sfiare the same equation and same 
initial data, but tie solutions differ as a consequence of a furtfier choice: of existence or 
not existence of a free boundary. 

• The pressure. There is a way of formally unifying the standard theory of free- 
boundary solutions for the porous medium equation with the present theory of free- 
boundaries in fast diffusion. This works by means of the new variable called pressure that 
for an equation of the form Vt = ^{v)xx is defined as 



If the integral is convergent at f = then the choice a = is made. Otherwise, any value 
in the domain of $ is good. Indeed, the similarity in the behaviour of the pressure of 
Porous Medium case, = u"^, m > 1, and our case (which lies in the power range 

m < 0) is seen when we write the partial differential equation for v which reads 



with a{v) = For = cm™ we get v = cmu"^ ^/{m — 1) and cr(i;) = {m — l)v. The 

paper |V03cj explores further the properties of free-boundary solutions for the pressure of 



jvEE|. 



(11.2) 




(11.3) 



Vt = a{v)Av + |Vf 
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fast diffusion equation. Note that, contrary to the porous medium case where the support 
of the pressure solutions expands in time, in our case it shrinks in time. 

• Grey levels as end values. Another variation of the main theme is to consider 
solutions losing the white or black level. This is investigated in the paper |RV95j . where 
many solutions of the differentiated equation for v are obtained in the range > m > — 1 
with the same initial data by assigning fluxes at infinity: 

(11.4) hm {v^-\,) = -fit), hm iv"'-'v^,) = g{t), 

X— >oo X— >— oo 

for bounded functions f,g>0. In terms of u it means in particular that we can obtain a 
free-boundary solution with a decreasing value of highest color level, 

m(oo, t) = F{t) = 1 - / fit) dt < 1. 

Jo 

Analogously, we can impose an increasing value for the lowest color level, 

u{-oo,t) = G{t) = [ g{t)dt>0. 
Jo 

• NoN-MONOTONE FRONTS. There is a gap in the theory we have developed, namely 
that we assumed that the front is monotone in the space direction, in other words, that 
u is nondecreasing as a function of x. The limitations of fast diffusion equations to admit 
non-monotone solutions are not accidental. In the paper |HVn2j it is proved that the 
model equation Ut = u^'^u^x does not admit non-monotone solutions if m < 0. In terms 
of the equation for v = u^ 'ii means that there are no solutions with both signs. However, 
there is hope when using functions $ like in ()1.4j) that are degenerate at infinity, i.e., 
$'(cxd) = 0, but regular at all other values, e.g., $ G C^(]R). Indeed, the problem with 
the monotonicity happens because of the singularity aX Ux = ^ which forbids maxima 
or minima, while our main interest is in the free-boundary that is governed by the large 
values, Ux oo. This question will be investigated separately. 

• NoN-MONOTONE NONLINEARITIES. Note that in our problem setting the equation 
is forward parabolic, and the backward movement of the interface is due to the effect 
of the singular boundary condition, which happens to be compatible with the equation 
for the appropriate class of functions $ (e.g., for powers m < 0). In the spirit of the 
Perona-Malik model, situations can be considered where $ is not monotone, hence the 
equation is not parabolic. A regularization is then needed so that one faces a regularized 
forward-backward diffusion problem. See in this respect |BBDUj and |Witj . This subject 
is again under study. 
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